Reduction of Pauli paramagnetic pair-breaking effect 
in antiferromagnetic superconductors 
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Antiferromagnetic superconductors in a magnetic field are studied. We examine a mechanism 
which significantly reduces the Pauli paramagnetic pair-breaking effect. The mechanism is realized 
even in the presence of the orbital pair-breaking effect. We illustrate it using a three-dimensional 
model with an intercalated magnetic subsystem. The upper critical field is calculated for various 
parameters. It is shown that the upper critical field can reach several times the pure Pauli param- 
agnetic limit. The possible relevance to the large upper critical field observed in the heavy fermion 
antiferromagnetic superconductor CePtaSi discovered recently is briefly discussed. We try to under- 
stand the large upper critical field in the compound CePtsSi and field-induced superconductivity in 
the compound CePb3 within a unified framework. 

PACS numbers: 74.25.-q, 74.25.Dw, 74.25.Op 



In this paper, we examine superconductivity in the 
presence of a coexisting antiferromagnetic long-range or- 
der. In a magnetic field, some types of antiferromagnet 
exhibit canted spin structures, which give rise to a net 
ferromagnetic moment proportional to the magnitude of 
the magnetic field. The ferromagnetic moment creates an 
exchange field on the conduction electrons by the Kondo 
interaction. When the Kondo interaction is antiferro- 
magnetic, the exchange field reduces the effective Zeeman 
energy. As a result, the Pauli paramagneticpair-breaking 
effect pi 0| can be significantly reduced 0, Q > an d the 
upper critical field H c2 can largely exceed the Pauli para- 
magnetic limit Hp (Chandrasekhar and Clogston limit). 

The analogous compensation mechanism in supercon- 
ductors with uniformly aligned spins is known as the 
Jaccarino-Peter mechanism 5] , which explains the field- 
induced superconductivity (FISC) in the compounds 
E ^ni^xMo 6 S 8 H, CePb 3 0, and A-(BETS) 2 FeCl 4 H 
l9l floL fll| ■ However, the resultant phase diagrams of the 
present mechanism are quite different from that of the 
Jaccarino-Peter mechanism. The superconductivity oc- 
cupies a large single area including the zero field (H = 0) 
in the T-H phase diagram, while depending on the pa- 
rameter values, the FISC may occur, where T and H 
denote the temperature and the magnetic field, respec- 
tively. In our previous papers 0, 3 , we have proposed 
this mechanism in a multilayer system in a parallel mag- 
netic field. We have adopted such a system, because the 
reduction mechanism is clearly illustrated. In this paper, 
we extend the same mechanism to more general systems 
in the presence of the orbital pair-breaking effect. The 
magnitude and the temperature dependence of the resul- 
tant upper critical field can be similar to those of the 
spin-triplet superconductivity of equal spin pairing. 

The present work was motivated by the heavy fermion 
antiferromagnetic superconductor CePtaSi discovered re- 
cently 12] . Due to the heavy quasi-particle mass, the 
orbital pair-breaking effect is much weaker than in con- 
ventional metals. Furthermore, the compound exhibits 



antiferromagnetic transition at Taf ~ 2.2 K and super- 
conducting transition at T c w 0.75 K at the zero field. 
The upper critical field H c2 (0) « 5 T at T = is much 
larger than the Pauli paramagnetic limit estimated by the 
simplified formula H P « 1.86[T/K] x T C (0) [K] « 1.4 T. 
Such a large H C 2 seems to suggest spin-triplet pairing, 
although the lack of inversion symmetry of this com- 
pound seems disadvantageous to spin-triplet pairing. Re- 
cently, Frigeri et al. have shown that spin-triplet pairing 
is not entirely excluded by the lack of inversion symme- 
try. Therefore, equal spin pairing may be the reason for 
a large H C 2 [13| . In this paper, however, we propose an- 
other scenario based on antiparallel spin pairing includ- 
ing spin-singlet pairing. Samokhin et al. have carried 
out band structure calculations and proved within their 
theory that the order parameter must be an odd func- 
tion of the momentum k [14j . We should note that this 
does not mean the occurrence of equal spin pairing [15| . 
Hence, our scenario does not contradict their result. 

In principle, we might be able to test these scenarios 
by rotating the magnetic field. The upper critical field 
H C 2 must be highly anisotropic, if the large H C 2 is due to 
equal spin pairing and the d-vector is fixed. For example, 
in the compound CePtsSi, the st rong spin-orbit interac- 
tion favors a particular d- vector [l3|. For the magnetic 
field parallel to the d-vector, the upper critical field must 
be strongly suppressed. In contrast, the upper critical 
field can be nearly isotropic in the present mechanism. 
The Knight shift measurement might also give useful in- 
formations |l5J . 

As mentioned above, the purpose of this study is to 
examine the present mechanism in the presence of the 
orbital pair-breaking effect. Therefore, we examine a 
three-dimensional system with an intercalated magnetic 
subsystem. We divide the magnetic subsystem into two 
sublattices, which we call A and B. We approximate the 
conduction electron system with an isotropic continuum 
system for simplicity. 

When we apply the theory to the compound CePtaSi, 



the A and B sublattices are alternate layers of Ce atoms, 
and the spins on each layer are ferromagnetically or- 
dered [ll| . The superconductivity and the antiferromag- 
netism may exist in the same degrees of freedom in this 
compound, but we simplify the situation by dividing the 
degrees of freedom into two coupled systems. We also ig- 
nore the spin-orbit interaction and the lack of inversion 
symmetry for simplicity, because they do not affect the 
present mechanism. 

First, let us describe the spin structure in the magnetic 
subsystem. We consider the situation in which the an- 
tiferromagnetic transition occurs at a temperature much 
higher than the superconducting transition temperature, 
which is consistent with the observation in the compound 
CePtsSi JI3. Since the ordered state is rigid, we neglect 
the modification of the magnetic structure by the occur- 
rence of superconductivity. Therefore, we consider the 
cases in which the magnetic subsystem is effectively de- 
scribed by the Hamiltonian 



Hs = X! "^ Si ' S J + ^ds^BH ■ Si 



(1) 



i,j 



where g s and fj,B denote the g- factor of the localized spins 
and the Bohr magneton |17|. We define z as the number 
of antiferromagnetic bonds for a given site. We set Jy = 
except for the nearest neighbor sites (i,j). For the 
compound CePtsSi, we assume that Jy = J11 < for 
(i,j) on the same layer (the same sublattice), while Jy = 
J± > for (i,j) on the adjacent layers (the different 
sublattices). In this case, we have z = 2. We could 
examine other crystal structures. For example, we have 
Z = 6 in the system on the cubic lattice. In any case, 
we assume that Jij = J > for the antiferromagnetic 
bonds. 

We take the x-axis of the spin space as being in the 
direction antiparallel to the external field. Hence, we set 
H = (—H,Q,0). The directions of the axes of the spin 
space do not necessarily coincide with the crystal axes 
in the present model. For T < T c <C Taf, we neglect 
the thermal fluctuations and regard the localized spins as 
classical variables. If we ignore the anisotropy energy of 
the localized spins in the magnetic layers, the sublattice 
magnetization appears in the direction perpendicular to 
the applied field, since it is energetically favored. Hence, 
we introduce the classical variable 9 by 



Si =S(sin0,O,±cos( 



(2) 



where we define the z-axis of the spin space as being in 
the direction of the sublattice magnetization. For the 
double sign ± in eq. 10, we adopt + and — signs when 
site i belongs to the A and B sublattices, respectively. 
Here, we have introduced the magnitude of the localized 
spins S taking into account the shrinkage of spins by 
fluctuations. The total energy of the magnetic subsystem 
is expressed as 



E(0) = NS 



1 JS cos 26* — 2h, sin( 



En 



(3) 



where Eq is a constant and h s : 
energy becomes minimum when 

zJSsmQ = 



g s H B H/2. The total 



(4) 



for h s /zJS < 1, i.e., H < 2zJS/g s /iB = Haf- For a 
high field in which H > Haf, the staggered moment dis- 
appears (0 = 7r/2). From now on, we mainly consider a 
weak field H < Haf, where the antiferromagnetic order 
persists. The localized spins are also described by the 
spin density 



S(r) = Y,5W(r-R i )S i 



(5) 



where Ri denotes the position of site i. 

We treat the degrees of freedom of the conduction elec- 
trons as a continuum model. The net spin moment S sin 6 
gives rise to the exchange field on the conduction elec- 
trons. In order to reproduce the exchange field, we em- 
ploy an extended Kondo Hamiltonian defined as 



H 



K 



d 3 r / dV J K (r-r') 



S(r)-J2bPUr')<T„*>^>(r% 



(6) 



where cr aa i and ipa{r) denote the Pauli matrices and 
the annihilation operator of the conduction electron at 
r. Here, we note that eq. 10 is used in the calcu- 
lation of the order parameter A(r), which appreciably 
varies in the length scale of the coherence length £ ^> a, 
where a denotes the lattice constant. Therefore, after 
inserting eq. JSJ in eq. JJjJ, we replace the spin vari- 
able Si by the spatial average (S)t t over the volume 
V£i on the order of £ 3 near site i, which is defined as 
( s )& = E«, e „ e< S J /Efl je „ 5 , 1- Thus > we obtain 



H K ^fd 3 r'J2JK(Ri-r') 



(^•EW( r V^( r ') 



(7) 



In actuality, the microscopic fluctuations omitted here 
may affect the superconductivity through the self-energy 
of the electrons. We regard those corrections to the nor- 
mal state as being included in the effective mass to*, if 
it exists. Since £ S> a, there are many localized spins in 
the volume v^i , the antiferromagnetic components of the 
spins vanish on average, while the ferromagnetic moment 
remains. Therefore, we obtain (S)^ « (S sin 9,0,0). For 
specific lattice structures, this equation is satisfied ex- 
actly. Furthermore, we assume that the localized spins 
on zk lattice sites around r' take part in the Kondo in- 
teractions with the conduction electron at r', and that all 
of these coupling constants are equal to Jr. Therefore, 
we obtain 



H 



K 



h, 



Jd 3 r'J2H(r> x ^^(r')], (8) 



where h ex denotes the exchange held h e 
From eq. (@J, we obtain 



g s /.i B z K J K 
2 zJ 



H, 



zkJkSsvuO. 



(9) 



for H < Haf- For H > Haf, we have h ex — zkJkS. 

Taking into account the exchange held h ex and the 
external held h, the kinetic energy term of the conduction 
electron Hamiltonian is written as 



H a 



E 



d 3 r V'i( r ) 



-h 2 



(l 



2m 



(V 



-A(r)) 2 S aa , 
c 



(10) 



where m* and h denote the band effective mass and the 
renormalized Zeeman energy h = h — h ex , respectively. 
Here, we have defined h = [i e H and [i e = g e /j,B/2 where 
g e denotes the g-factor of the conduction electrons. For 
H < Haf, we have 



~ h=(l _£sZKJK )h= ( 1 _ JK 



g e zJ 



J 



(11) 



.4/-' 



obtained by eq. 
J'af = g e zJ/g s z K - 



(J5J), where we have dehned 
In contrast, for H > Haf, we 
have h = h — zkJkS 

In eq. (|ll|l . it is explicit that the Zeeman energy is 
reduced by the exchange held when Jk > 0. It is found 
that the reduction effect becomes maximum when Jk — 
J'af- Obviously, this reduction mechanism is not affected 
by the anisotropy of the order parameter. Hence, we 
consider the s-wave superconductor as an example. 

Applying the standard theory of superconductivity to 
the present system, we obtain the integral equation for 
the superconducting transition temperature T c |l8l Il9l 



In- 



T r 



n 



(0) 



i r°° i r 

- \ dx / d9sm9 

2 J smh x J 



exp( x 2 sin 2 6*) cos(2h'x) — 1 



(12) 



,(0) 



where T c denotes the zero-field transition tempera- 
ture, and we have defined k! = (vf /2irT c ) 2 2\e\H / 'c and 
h! = h/2nT c . Here, we have ignored the possibility of 
the Fulde-Ferrell-Larkin-Ovchinnikov state. The exten- 
sion to include it is straightforward [j, |2j|. From the 
form of k! . it is convenient to define a constant 



v F 



'2ttT, 



(0) 



)2 2H 
c 



(13) 



with which we can write k' = a m H(Tc /T c ) . We also 
define the ratio of the scales of the paramagnetic and 
orbital effects as 



h/2irT ( 
a m H 



(0) 



g e irTc m* 
8cf rn 



(14) 
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FIG. 1: Temperature dependence of upper critical fields for 
r m — 0.8. The dotted curve shows the pure Pauli param- 
agnetic limit |23|| in the absence of the orbital pair-breaking 
effect at Jk = 0, where g e = 2. The curves of the upper crit- 
ical field for Jk/J'af = IT, 1.2, and 1.5 coincide with those 
for Jk/J'af = 0-9, 0.8, and 0.5, respectively. 
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FIG. 2: Critical fields of superconductivity for Jk = 1.2J' AF , 
r m — 0.8 and a m H cent = 1. The dotted and short dashed 
curves show the pure paramagnetic and pure orbital limits, 
respectively. The dot-dashed line shows the critical field of 
the antiferromagnetic phase. 



where ep = m*v 



,/2- 



Equation (|12|) is solved numerically. Figure^shows the 
results for r m = 0.8 and Jk < J'af- H i s found that the 
upper critical field is enhanced as Jk increases up to J' AF ■ 
When Jk i= 0, the upper critical field H c2 can exceed 
the Pauli paramagnetic limit Hp(T). In particular, for 
Jk ~ J'af i ^ ne u PP er critical field H cZ can reach fourfold 
the Pauli limit Hp. When Jk = J'af, t ne pure orbital 
limit H C 2o(T) is recovered. For larger r m , i.e., weaker 



orbital effect, the ratio H C 2/Hp increases. 

If there is a temperature region in which H c2 (T) > 
Haf in Fig. ^ the curve of H C 2(T) should be modified 
there so that it saturates more rapidly, because the ex- 
change field is constant for H > Haf- In the scale of 
the right vertical axis, H = Haf gives r m a m HAF — 
g e zJS/2-Kg s Tf\ Hence, if g e /g s » 1 or if zJ5/T c (0) ~ 
T AF /Tf> > 1, the upper critical field H c2 {T) does not 
exceed Haf in the whole temperature region. 

In a high-field region in which H > Haf, where the 
spins are aligned uniformly, the mechanism is reduced to 
the Jaccarino-Peter mechanism. In this case, the Zeeman 
energy completely vanishes at H = zkJkS ' / [i e = H cent . 
Therefore, if H cent + Hp > Haf is satisfied, we need 
to examine the possibility of FISC. Since H^t/HAF = 
Jk/J'afi FISC occurs only when J K > J' AF 3]. Figure^ 
shows the critical fields for Jk — ^-2J'af-> r ™, = 0.8 and 
a m H cent — I. This set of the parameter values gives 
a m HAF ~ 0.83, and if g s = g e = 2, Taf ~ zJS — 
2irg s Tc r m a m HAF 1 9e ~ 4T^ . In Fig. [2J we find an 
area of FISC inside the region of H c20 (T) > H > Haf- 

In this context, the compound CePb3 corresponds 
to the case of H cent ~3> Haf, i-e., Jk ^ J'afi an d 
r m ^> 1 |2J|. In contrast, the compound CePt 3 Si has 
two possibilities, i.e., Jk ^ J'af an< ^ J~K ~ J'af- For 
the former case, FISC might be observed at very high 
fields, if H c20 (0) > H cent > H AF .25]. In the com- 
pound UPd2Al 3 , it was observed that Taf — 14.3 K 
and T c (0) = 2.0 K |2Jj. Such a high T AF suggests small 
Jk/ J'af-i so that the reduction of the paramagnetic effect 
is not pronounced. This coincides with the experimental 
fact that H C 2 < Hp and the absence of FISC in the com- 
pound UPd2Al3. Generally speaking, the difference in 
the experimental results Taf ~ 1.1 K, 2.2K, and 14. 3K 
in CePb 3 , CePt 3 Si, and UPd2Al 3 , respectively, is consis- 
tent in the present theory with their phase diagrams, if 
their Jk 's are of the same order. 



Here, the condition r m ~ 1 can be satisfied in heavy- 
fcrmion superconductors, but not in conventional metals, 
because r m ex g e Tc /wlvf 2 — g e {Tc /eF)(m*/m). Also 
in the scenario with equal spin pairing, we need r m of 
order 1, in order to reproduce a large ratio of H c2 / Hp 
as observed. 

In conclusion, we have shown that the Pauli param- 
agnetic pair-breaking effect can be considerably reduced 
in antiferromagnetic superconductors, even in the pres- 
ence of the orbital pair- breaking effect, unless r m <C 1. 
One of the necessary conditions for the occurrence of the 
present mechansim is that the Kondo exchange coupling 
constant Jk is positive and comparable to the scale of 
the antiferromagnetic exchange energy J'af- The present 
result may explain the high H c2 observed in the com- 
pound CePt 3 Si. The phase diagrams of the compounds 
CePt 3 Si and CePb 3 can be understood within a unified 
framework. The positive Jk can be attributed to su- 
perexchange or kinetic exchange processes. Hence, the 
condition Jk ~ J ~ J'af can be satisfied, because the 
antiferromagnetic exchange interaction J originates from 
exchange processes similar to those for Jk in the same 
crystal structure. The hybrid ruthenate-cuprate com- 
pound RuSr 2 GdCu 2 8 HJ [H, |5jj is another possible 
candidate of the present mechanism, if one applies a par- 
allel magnetic field to it. The antiferromagnetic long- 
range order with weak ferromagnetism due to the canted 
spin structure has been observed in this compound at 
zero field |I|. 
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